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any known SGR. Subsequent observations of the source error box with the 
Swift/X-Taj Telescope and the Rossi X-ray Timing Explorer {RXTE) led to 
the discovery of a new X-ray source, with a spin period of 7.56 s, confirming 
SGR J1833— 0832 as a new magnetar. Based on our detailed temporal and spec- 
tral analyses, we show that the new SGR is rapidly spinning down (4xl0~^^ 
s/s) and find an inferred dipole magnetic field of 1.8x10^^ G. We also show that 
the X-ray flux of SGR J1833— 0832 remained constant for approximately 20 days 
following the burst and then started to decline. We derived an accurate loca- 
tion of the source with the Chandra X-ray Observatory and we searched for a 
counterpart in deep optical and infrared observations of SGR J1833— 0832, and 
for radio pulsed emission with the Westerbork Radio Synthesis Telescope. Fi- 
nally, we compare the spectral and temporal properties of the source to the other 
magnetar candidates. 

Subject headings: pulsars: individual (SGR J1833— 0832 ) — X-rays: bursts 



Introduction 



The first c onfirmed Soft Gamma Repeater (SGR) burst was detected in 1979 but it 
took ~ 7 years (JLaros et al.lll987l : lAtteia et al.lll987l : iKouveliotou et al.lll987l ) to associate it 
with a newly discovered population of astronomical objects, namely neutron stars repeatedly 
emitting soft and short gamma-ray bursts with peculiar properties. A decade later, these 
objects, together with their cousins. Anomalous X-ray Pulsars (AXPs; see e.g., iMereghetti 
20081 ). were shown to possess extremely high magnetic fields (~ 10^^ G), as deduced from 
their spin periods and spindown rates, and were named magnetars (iDuncan fc Thompson 
I992I : [Thompson fc Dunca3ll996l : IKouveliotou et"aDll998l : IVasisht fc Gotthelf Ill997l ). The 
most distinguishing magnetar characteristics are a narrow range of spin periods (2— 12s), 
high spindown rates (~ 10~^° s/s) and unpredictable burst-active episodes. The latter 
can have a multitude of bursts with a wide luminosity range, anywhere from 10~^ to 10^ 
times the Eddington luminosity. It is these bursts that enable the initial classification of a 
source as an SGR, to be subs equently confirme d as a magnetar with timing analysis (see 
Woods fc Thompson! I2OO6I and lMereghettill2008l for a detailed review on magnetars.) 



By 2008 a dozen sources were confirmed magnetars and a couple more were listed as 
magnetar candidates pending observations of new burst active episodes. For the known 
sources, frequent spin period monitoring observations with the Rossi X-ray Timing Explorer 
[RXTE) over the last fifteen years have provided a wealth of insight about their dynamics and 
energetics. Interestingly, despite their common characteristics, there are several important 
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differences across this small population, and a larger sample of objects at different evolution- 
ary stages is necessary to better understand the phenomenon. Magnetar hunting however, 
is not an easy task; we have two strikes against us ab initio: their very small numbers and 
their very short active periods (days to months). Thus, the best way to discover magnetars 
is the combination of an all-sky X-and gamma-ray monitor - to detect their bursts - and a 
timing instrument - to measure their spin and spindown rates, and estimate their magnetic 
fields. 

On 2010 March 19, the Burst Alert Telescope (BAT) onboard NASA's Swift satellite 
triggered at 18:34:50 UT on one short and soft burst (Barthelmy et al. 2010) originating from 
a region in the plane of our Milky Way galaxy that harbors at least three other magnetar 
candidates: AXP XTE J1810-197, and SGRs 1806-20 and 1801-23. Using the initial source 
localization with the BAT, which had a 3' error radius, we were able to excluded all these 
sources; further, based on the spectral and temporal properties of the event, we identified the 
source of the emission as a possible new SGR. Subsequent observations (starting from 60 s 
after the BAT trigger) of the source region with the Swift /X-Ray Telescope (XRT) provided 
the source position with a 1.6" uncertainty ( iGelbord et al.l 120101 ). 5'w/i/Ultraviolet Optical 
Telescope (UVOT) observations starting 71 s afte r the trigger did not detec t a source in 
the UV or optical bandpass at the XRT position ([Marshall &: Gelbordll2010l ). To confirm 
the nature of the source, we observed it with RXTE for 1.3 ks less than 3.25 hours after its 
discovery was reported. A subsequent timing analysis of this ob servation revealed c oherent 
pulsations in the X-ray flux of the source with a period of 7.56 s (JGogiis et al.ll2010al ). Thus, 
we were able to identify SGR J1833— 0832 as a slowly rotating neutron star in less than 
6 hours after its discovery, the second time such a fast discovery was enabled only due 
to the unpre cedented synergy o f Swift and RXTE (the first was the identification of SGR 
J0501+4516; iGoeiis et allboooh . 



Swift was the only satellite that triggered on SGR J1833— 0832 ; although the source was 
in the field of view of the Gamma ray Burst Monitor (GBM) onboard the Fermi Gamma-ray 
Space Telescope (FGST), the fluence of the event was much lower than any of the GBM 
trigger thresholds. To precisely locate the source we triggered our Target of Opportunity 
(ToO) observations with the Chandra X-Ray Observatory [CXO) and obtained observations 
with the European Southern Observatory (ESQ) Very Large Telescope/High Accuity Wide- 
field K-band Imager ( VLT/HAWKT) and with the Gran Telescopio Canarias/ Optical Imager 
and Spectrograph (GTC/OSIRIS). We also observed with the Westerbork Synthesis Radio 
Telescope Pulsar Machine II ( WSRT/PuMa IT) and searched for radio pulsations from the 



source. 



In this paper, we report on the discovery of SGR J1833— 0832, and introduce its char- 
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acteristics through our detailed multi-wavelenth investigations. In §2, we describe all multi- 
wavelength data used in our study. We present in §3 the precise localization of the new 
source using CXO, 2MASS and UKIRT data, and the temporal and spectral characteristics 
of the persistent X-ray emission from the source using our Swift and RXTE observations. 
Finally, we discuss and interpret our results in §4. 



2. Observations 
2.1. Swift/BAT 



Fo llowing the 2010 March 19 BAT discovery of a burst (trigger=416485; iGelbord et al. 



(I2OIOI ): iBarthelmy et al.l (120101 ). Swift slewed immediately to the on-board BAT position of 
RA, Dec = 18'^33'"46^ -08°32'13" (/, h = 23^^19"^38^ 0° 0' 22"; radius uncertainty of 3', 90% 
containment). This location is 13.3° from the position of SGR 1806—20. The Swift/XBT 
detected a new X-ray source within the BAT error circle; these observations are described 
in §2.21 This was the only triggered burst detected from the source with BAT and we report 
below on its spectral and temporal characteristics. 

The burst was detected within the 87% BAT partial coding region; its mask-weighted 



light curve (Figure [T]) shows a single peak structure. The T90 duration (JKouveliotou et al. 



I993I ) of the burst is 0.013 ± 0.002 s (15 — 350 keV; total error including systematics) . 



We found that several continuum models could provide a statistically acceptable fit to 
the time-integrated spectrum (Tq to To+0.016 s); a power law with an exponential cutoff u 
with a photon index, a, of —1.02 ±1.35 and an -Epcak of 38.2 ±5.0 keV (x^/degrees of freedom 
(d.o.f.) = 43.5/56), a blackbody with kTBB = 10.0 ± 1.2 keV (xVd-o.f. = 43.8/57), and an 
optically thin thermal bremsstrahlungj with /cTottb = 29±5 keV (x^/d.o.f = 52.2/57). For 
the cutoff power law model, the total fiuence in the 15 — 150 keV band is (1.3 ± 0.2) x 10^^ 
erg/cm^ and the 6 ms peak flux (measured from Tq ± 2 to Tq ± 8 ms) is (9.9±2.1)xl0^^ 
erg/cm^/s (19.4 ± 3.8 photons/cm^ /s). All quoted errors are at the 90% confldence level. 



^Power law with an exponential cutoff is parametrized as dN/dE 
A{E/hQ keV)-" exp {-E{2 ~ a)/Speak) 

^OTTB is parametrized as dN/dE = ^(i^/SOkeV)-^ exp {~E/kT) 
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Fig. 1.— The Swift/BAT lightcurve of the burst from SGR J1833-0832 in the 15 - 350 keV 
band. The time bin size is 2 ms. 



2.2. Swift/XRT 

Swift/XBT observations of SGR J1833— 0832 were obtained as a ToO monitoring pro- 
gram starting on 2010 March 19 at ~ 10 ks per observation for a total exposure of ~184 
ks (Table [1]). The XRT observations were coordinated with RXTE in a manner that com- 
plemented and optimized usage of both instruments. The XRT data were accumulated in 
Windowed Timing (WT) mode only during the third observation (002, see Tabled]), while 
all other data were collected in Photon Counting (PC) mode. In the WT mode only the 
central 8' of the field of view is read out, providing one-dimensional imaging and full spectral 
capability with a time resolution of 1.8 ms. The PC mode provides full spatial and spectral 
resolution with a time resolution of 2.5 s. 

The data were processed with standard procedures (xrtpipeline vO.12.1) and filtering 
and screening criteria, using FTOOLS in the Heasoft package (v. 6. 6.1). We used a 0-2 and 
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0-12 grade selection for data in the WT and PC modes, respectively. This selection provides 
the best combination of spectral resolution and detection efficiency. SGR J1833— 0832 was 
imaged far from the CCD hot columns so that no corrections inside the photon extraction 
regions were necessary. 

We extracted the source events for timing and spectral analysis from a circular region 
of 47" radius (20 pixels) centered on the source position as determined with xrtcentroid. 
The background was extracted from an annular region centered on the source with 90" 
and 120" internal and external radii, respectively. Ancillary response files were generated 
with XRTMKARF, and account for different extraction regions, vignetting, and Point Spread 
Function corrections. We used the spectral redistribution matrices vOll in CALDB. All 
spectral fits were performed using XSPEC v. 12. 5. 10 (Arnaud 1996). 



2.3. RXTE 

We monitored the source with the i?Xri?/Proportional Counter Array (PCA) starting 
on 2010 March 19, 21:49:10 UT (i.e., ~3.25 hrs after the BAT trigger) until 2010 May 6 in 
31 pointings. The time spacing between consecutive observations was on average 1.6 days, 
ranging between 0.06 and 4.9 days. The exposure of individual observations varied between 
1.3 ks and 12 ks, with a total exposure of 178.6 ks over 47 days. In Table [21 we list a log of 
RXTE observations of SGR J1833-0832. 

In this study we employed PCA data collected in the 2 — 10 keV range. We inspected the 
lightcurve of each observation binned with 0.03125 ms to search for additional SGR bursts, 
and identified four events with duration <12 ms. These events were rather weak compared 
to the typical bursts seen with the PCA from other SGRs; their detection indicates that a 
low level of bursting activity continued after the BAT trigger. We then filtered out the times 
of these four short bursts to obtain burst-free event lists for our timing studies. Finally, we 
converted all event arrival times to the Solar System Barycenter. 



2.4. Fermz/GBM 

Although GBM did not trigger on the Swift/BAT burst, we searched the daily GBM 
data (with time resolution of 256 ms in continuous naode a nd 64 ms in trigger mode) using 



our untriggered event search algorithm (JKaneko et al.ll2010l ) to uncover any additional weak 



events from the source. The search was performed in the 10—300 keV band on a 20-day long 
segment starting from 2010 March 15 (i.e., 4 days prior to the BAT trigger). We did not 
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Table 1: Summary of the Swift/XBT observations. 



Seq. 


Start time (UT) 


End time (UT) 


Exposure 

(s) 


000 


2010-03-19 18:35:50 


2010-03-20 15:38:07 


29174 


001 


2010-03-21 04:14:14 


2010-03-21 15:41:31 


10756 


002 


2010-03-21 16:47:37 


2010-03-22 02:32:06 


9954 


003 


2010-03-22 07:50:12 


2010-03-22 20:48:15 


9869 


004 


2010-03-23 08:58:47 


2010-03-23 23:39:46 


8068 


005 


2010-03-24 01:07:40 


2010-03-25 00:00:13 


10266 


006 


2010-03-25 01:36:34 


2010-03-26 00:00:12 


9844 


007 


2010-03-26 12:28:42 


2010-03-26 22:33:11 


9961 


008 


2010-03-27 10:57:33 


2010-03-27 21:02:31 


9855 


009 


2010-03-28 00:04:08 


2010-03-28 18:02:39 


10919 


010 


2010-03-29 12:44:13 


2010-03-29 21:13:35 


9511 


Oil 


2010-03-31 00:20:50 


2010-03-31 23:11:58 


7942 


012 


2010-04-03 14:59:40 


2010-04-03 23:29:26 


10028 


013 


2010-04-07 13:45:23 


2010-04-07 22:15:53 


10097 


014 


2010-04-12 09:41:24 


2010-04-12 22:42:57 


5130 


015 


2010-04-13 03:08:11 


2010-04-13 08:21:51 


4054 


016 


2010-04-15 00:01:07 


2010-04-15 11:23:58 


9417 


017 


2010-04-18 14:53:02 


2010-04-19 07:10:56 


8817 



Table 2: Summary of the RXTE observations. 



ObsID 


Start Time (UT) 


End Time (UT) 


Exposure (s) 


95048-03-01-10 


2010 Mar 19 21:49:10 


2010 Mar 19 22:19:10 


1335 


95048-03-01-00 


2010 Mar 21 08:07:28 


2010 Mar 21 12:13:14 


10443 


95048-03-01-01 


2010 Mar 21 12:50:24 


2010 Mar 21 15:22:10 


6892 


95048-03-01-02 


2010 Mar 21 17:32:32 


2010 Mar 21 23:13:10 


13396 


95048-03-01-04 


2010 Mar 23 05:56:32 


2010 Mar 23 11:03:10 


11936 


95048-03-01-03 


2010 Mar 23 15:03:28 


2010 Mar 23 19:08:10 


9719 


95048-03-01-05 


2010 Mar 24 08:19:28 


2010 Mar 24 09:17:10 


3425 


95048-03-01-07 


2010 Mar 24 11:27:28 


2010 Mar 24 12:19:10 


3071 


95048-03-01-08 


2010 Mar 25 07:52:32 


2010 Mar 25 10:21:10 


6723 


95048-03-02-00 


2010 Mar 26 04:25:20 


2010 Mar 26 06:49:10 


6451 


95048-03-02-01 


2010 Mar 27 03:55:28 


2010 Mar 27 06:21:10 


6562 


95048-03-02-02 


2010 Mar 28 15:59:28 


2010 Mar 28 18:18:08 


5412 


95048-03-02-03 


2010 Mar 29 23:19:28 


2010 Mar 30 01:39:10 


5960 


95048-03-02-04 


2010 Mar 30 22:51:28 


2010 Mar 31 01:02:10 


5616 


95048-03-02-05 


2010 


Apr 1 06:27:28 


2010 Apr 1 07:12:10 


2633 


95048-03-02-06 


2010 


Apr 1 07:50:24 


2010 Apr 1 08:45:10 


3283 


95048-03-03-00 


2010 


Apr 2 02:39:28 


2010 Apr 2 05:10:08 


6804 


95048-03-03-01 


2010 


Apr 3 00:37:20 


2010 Apr 3 01:34:10 


3341 


95048-03-03-06 


2010 


Apr 3 02:11:28 


2010 Apr 3 03:08:10 


3371 


95048-03-03-02 


2010 


Apr 3 21:01:20 


2010 Apr 3 23:32:10 


6774 


95048-03-03-05 


2010 


Apr 6 19:38:24 


2010 Apr 6 20:35:10 


3321 


95048-03-04-00 


2010 


Apr 9 02:33:20 


2010 Apr 9 04:34:48 


6810 


95048-03-04-01 


2010 


Apr 13 00:58:24 


2010 Apr 13 03:01:10 


5116 


95048-03-05-00 


2010 


Apr 18 21:55:28 


2010 Apr 19 00:26:10 


6824 


95048-03-06-00 


2010 


Apr 23 18:03:28 


2010 Apr 23 20:35:10 


6949 


95048-03-06-01 


2010 


Apr 26 15:05:20 


2010 Apr 26 15:53:10 


2831 


95048-03-06-02 


2010 


Apr 26 16:40:32 


2010 Apr 26 17:36:10 


3309 


95048-03-07-00 


2010 


Apr 30 13:14:24 


2010 Apr 30 15:45:10 


6892 


95048-03-07-03 


2010 


May 3 11:50:24 


2010 May 3 12:46:08 


3279 


95048-03-07-01 


2010 


May 3 13:24:32 


2010 May 3 14:21:10 


3347 


95048-03-07-02 


2010 


May 6 10:27:28 


2010 May 6 12:58:10 


6773 
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find any untriggered events originating from the direction of SGR J1833— 0832, including 
the BAT-triggered burst and the bursts detected with RXTE. This is most hkely due to the 
extremely short duration of the bursts (< 16 ms), which made them practically undetectable 
over the background with the coarse time resolution of the GBM continuous data. We also 
inspected the GBM data for the times of the four short and weak RXTE bursts and again 
found no evidence of activity. 



2.5. Chandra 

We initiated ToO observations with the Chandra Advanced CCD Imaging Spectrometer 
(ACIS-I) in the Timed Exposure (TE) mode on 2010 March 23 for an effective exposure 
of 33 ks. The source was relatively bright at the time of the ACIS-I observation, resulting 
in a (background-subtracted) count rate of 0.104 =b 0.003 counts/s (0.5 — 8 keV), which 
is expected to cause a pile-up at the ~13% level. Our instrument choice was based on 
our goal, which was primarily to determine an accurate X-ray position of the source, as 
described in §3.11 We selected data collected in the 0.5—8.0 keV band with all four ACIS-I 
chips and generated a binned image (by a factor of 2) to improve the signal to noise ratio. 
We then employed wavdetecto to search for sources with a detection threshold of 5a above 
background and identified six uncatalogued point sources, including SGR J1833— 0832. We 
then repeated the same source finding procedure with an unbinned image to determine the 
accurate X-ray position of the SGR, which was at the aim point of ACIS-I. We list the X-ray 
positions ( J2000) of all sources in Table |3l The last column in the table indicates whether 
an IR counterpart was identified for a source in the Two Micron All Sky Survey (2MASS) 
data (see also §3.ip . 



2.6. Optical and IR Observations 

The field of SGR J1833— 0832 was observed with the imager OSIRIS, mounted on the 
10.4m Gran Telescopio Canarias (GTC) telescope, starting on 2010 March 20 at 05:41:13 
UT (i.e., 11.11 hrs after the BAT trigger). The observations were performed under non- 
optimal atmospheric conditions, with an average seeing of ~ 1.6", and consisted of 20x90 
s long exposures of a 1.85' x 2.35' portion of the sky using a Sloan z-band filter. The left 
image of Figure |2] shows the combined GTC image of all 20 exposures of the field around 
SGR J1833-0832. The mean time of the observations is 06:04:04 UT on 2010 March 20. 



^A source finding algorithm in CIAO 4.0 
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Table 3: X-ray point sources in the Chandra ACIS— I field of SGR J1833— 0832. 



Source 


RA (J2000) 


Dec (J2000) 


S/N 


2MASS 


Name 


(hh mm ss.ss) 


(dd mm ss.ss) 


Ratio 


Association 


CXOU J183344.4-083108 


18 33 44.38 


-08 31 07.71 


53.7 


SGR, No 


CXOU J183320.2-083426 


18 33 20.17 


-08 34 26.34 


6.1 


Yes 


CXOU J183322. 7-082409 


18 33 22.72 


-08 24 08.89 


6.5 


Yes 


CXOU J183317.9-082516 


18 33 17.86 


-08 25 15.45 


6.1 


Yes 


CXOU J183316.2-083615 


18 33 16.22 


-08 36 15.22 


8.7 


Yes 


CXOU J183359.4-082229 


18 33 59.44 


-08 22 29.18 


6.4 


No 



We do not detect any optical source within the Chandra error circle in the combined image 
down to a 3a limiting magnitude of 24.9. 

We also observed the source using the ESQ VLT/B.AWKI in the Ks-band under excellent 
seeing conditions (0.45"). A total of 44 dithered exposures of duration 60 seconds were 
obtained, providing a total exposure time of 2640 s. The mean time of the observations was 
09:01:09 UT on 2010 March 25. There is no obvious IR source within the Chandra error 
circle; further astrometry and limiting magnitude information is given below ( §3.11) . 



2.7. WSRT/PuMa II Radio Observations 

We observed SGR J1833— 0832 on three epochs in March, all within a week of the source 
discovery, using WSRT at central frequencies of 2280 and 1380 MHz (see Table H] for details). 
The pointing position was RA = 18'^33'^44^4, Dec = -08°31'08".4. WSRT was used in its 
tied-array mode, but as the position of the source was relatively uncertain, only telescopes 
out to baselines of 1.5 km were used. Furthermore, since one telescope was unavailable, 
only 11 telescopes were combined, giving a resolution of approximately 15" and 24" for 
2280 and 1380 MHz, respectively, and a gain of ~1 K/Jy. The WSRT splits the 160 MHz 
bandwidth into 8 bands of 20 MHz wide, and each band is recorded at baseband using the 



PuMa //pulsar backend ( jKaruppusamy et al.ll2008l ). The baseband data were converted into 



filterbank format with 64 channels per 20 MHz band and a sampling time bin of 8.192 ms. 
The resulting 512-channel filterbank data were searched blindly for a periodicity matching the 
known period of SGR J1833— 0832, and also folded using an ephemeris derived from the X-ray 
data ( §3.2p . We searched over a range of dispersion measures of 0-1800 pccm~^, where the 
maximum is that expe cted along the line of sight through the Galaxy in the NE 2001 model 



(jCordes fc Lazidl2002l ). The long spin period of the source means that special attention had 
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to be given to identify potential radio frequency interference. However, while there were some 
long period features present in the data, they did not significantly affect our sensitivity to 
the source. Inspection of both the folded pulse profiles and the candidate lists from the blind 
searches did not reveal any s ignificant detections of SGR J1833— 0832. Using the radiometer 



equation (JDewey et al.lll985l ) and considering pulse widths of 50% and 10%, which span the 
range of likely pulse widths for radio emitting magnetars, we obtain sensitivity limits of 0.6 
and 0.2 mJy at 2280 MHz and 0.3 and 0.1 mJy at 1380 MHz, corresponding to the 50% and 
10% widths, respectively. 

Table 4: WSRT - PuMa II Observation Log 



Epoch 


MJD (day) 


Duration 


(hour) 


Frequency (MHz) 


23 March 2010 


55278.1570 


1 




2280 


23 March 2010 


55278.1985 


1 




1380 


24 March 2010 


55279.1962 


1 




1380 


26 March 2010 


55281.2034 


1 




1380 



3. Data Analyses and Results 
3.1. Precise Source Location and IR Counterpart upper limits 

3.1.1. Imaging and astrometry 

We searched the 2MASS archival data for IR counterparts of the five X-ray sources 
significantly detected in the Chandra image ( §2.5p . in addition to SGR J1833— 0832 itself. To 
perform astrometry, we utilize observations taken as part of the United Kingdom Infra- Red 
Telescope (UKIRT) Deep Sky Survey (UKIDSS), which obtains deeper images than 2MASS 
(limiting magnitudes of K~ 18) with an improved PSF to aid with the elimination of source 
confusion. We show this improvement graphically in Figure [3l where we plot 2MASS and 
UKIRT images of each of the Chandra X-ray sources (excluding SGR 1833-0832, which is 
shown in Figure [2]). We clearly see that one source, confused in the 2MASS images (CXOU 
J183320.2-083426), can be resolved by UKIRT; a second source (CXOU J183317.9-082516) 
remains apparently confused in the UKIRT images due to the contribution of a bright and 
faint source within the Chandra position. For the two remaining sources (CXOU J183320.2- 
083426 and CXOU J183316.2-083615) we confirm associations with 2MASS sources. To 
obtain astrometry of the field we first aligned the UKIRT data to 2MASS, which gives a fit 
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with an RMS of ~ 0.20". Because of the source confusion in the 2MASS images at faint 
magnitudes, we hmit ourselves to stars in the magnitude range 8 < -R' < 12. We then 
compared the predicted locations of the X-ray sources with their measured centroids, and 
fit the resulting data with a simple linear shift in x and y coordinates (since we do not 
have a large number of sources for comparison this method provides a better estimate than 
allowing a larger number of free parameters). The resulting shifts are small: -0.25" in RA 
and 0.2" in Dec. The position of SGR J1833-0832, accurate to 0.40" is, RA = 18^33"^44^37, 
Dec = -08°31'07".5 (J2000). 




Fig. 2.— A zoomed in region around SGR J1833-0832 in the GTC/OSIRIS image (left) and 
in the VLT/HAWKI image (right), showing the C/ianc^ra position (Icr error radius). Labels 
A and B indicate the objects close to the Chandra position of SGR J1833— 0832 (see §3.1.2). 



This accuracy, however, is dominated by the error in alignment of the UKIRT im- 
ages to an absolute World Coordinate System (WCS; in this case 2MASS). A more accu- 
rate result was obtained by performing simple relative astrometry between the UKIRT and 
the Chandra imaging, providing a position accurate to ~ 0.3" (using the relative location 
of SGR J1833— 0832 on the UKIRT frames). Although there is no obvious counterpart 
for SGR J1833— 0832 in the UKIRT images, our deepest plausible limits come from our 
VLT/HAWKI observation with a total exposure of 2640 seconds. The very small field of 
view (FoV) of these observations precluded performing astrometry with Chandra sources di- 
rectly (none of the five X-ray sources were in the HAWKI FoV). However, we could perform 
relative astrometry between the UKIRT and the VLT images with a minimal error (<0.05"), 
and hence place SGR J1833-0832 on the VLT/HAWKI images with ~ 0.30" accuracy. The 
resulting location of SGR J1833— 0832 is shown in Figure O 
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J1 83320.2-083426 J1 83322.7-082409 J1 8331 7.9-08251 6 J1 8331 6.2-08361 5 J1 83359.4-082229 

Fig. 3. — Chandra locations for five X-ray sources overlayed on UKIRT images (top) and 
2MASS images (buttom). As can be seen the UKIRT allows an additional two counterparts 
to be identified (CXOU J183320.2-083426 and CXOU J183322.7-082409), which were am- 
biguous in the 2MASS data. 

3.1.2. IE limits 

We do not detect a source within the la Chandra error circle of SGR J1833— 0832 in 
the VLT /HAWKI m\a,ges. Performing photometric calibration relative to 2MASS yields an 
upper limit for the SGR J1833— 0832 counterpart within the C/iandra/ ACIS-I error circle as 
Ks > 22.4 {3a). We should note, however, that this is an idealized scenario since confusion 
between faint sources is an issue at the lowest flux levels. While there are no sources within 
the la Chandra error radius, there are two sources at the ~ l.Scr level from the Chandra 
position of SGR J1833— 0832: a moderately bright source (marked A in Figure [2]) with 
Ks = 18.2 ± 0.1, and a fainter source (marked B) with Kg = 21.2 ± 0.2. Since neither fall 
within the Chandra location, and due to the crowded nature of the field, we cannot make 
any strong claims for their association with SGR J1833— 0832. Future observations to search 
for variability of these sources may provide a stronger diagnostic. 



3.2. Pulse Timing Analysis 

We used an epoch folding algorithm to determine the pulse ephemeris of SGR J1833— 0832, 
as follows. We first determined the source pulse frequency using closely spaced RXTE/PCA 
observations on 2010 March 21 with a total exposure of 44 ks and generated a template 
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pulse profile in the 2—10 keV band. We then grouped all available PCA observations into 22 
discrete time segments, each separated by at least 0.3 days. For each segment, we generated 
the pulse profile and cross-correlated it with the template to obtain the phase drift with re- 
spect to the template. We finally fit the phase shifts with a first or higher order polynomial 
to determin e the pulse ephemer is of the source. The method summarized here is described 



in detail by IWoods et al. 



(120021 ). 



We first attempted to model the phase drifts with respect to the template using a 
first order polynomial. This resulted in a poor fit (x^/d.o.f. = 122.1/21) with significant 
deviations of phase drift values from the model. Modelling the phase drifts with a second 
order polynomial results in a good fit (x^/d.o.f. = 19.9/20) with a spin period of 7.5654091(8) 
s and a period derivative of (4.39±0.43)xl0~^^ s/s (epoch MJD 55276.5). We present in 
FigureHlthe pulse phase drifts with respect to the phase of the template profile with a second 
order polynomial fit (top panel), and the residuals of the fit (lower panel). 



3.3. Pulsed Fraction and Pulsed Intensity History 

We investigated the energy dependence of the RMS pulsed fractioiuof SGR J1833— 0832 by 
folding all available XRT data in PC mode in the 0.2—4 keV, 4—6 keV and 6—10 keV bands 
using the pulse ephemeris presented in §3.2. In Figure 0, we present the pulse profiles of 
SGR J1833— 0832 in the three specified energy intervals (left) and the normalized profiles 
obtained by dividing with the maximum value of each pulse profile (right). We note that 
although the intensity of pulsations varies significantly in energy, the normalized amplitude 
of the pulsed signal remains remarkably constant. Consequently, we find that the RMS 
pulsed fraction of SGR J1833— 0832 does not vary as a function of energy with RMS pulse 
fractions of (34±2)%, (36±2)% and (39±5)% in the above energy intervals, respectively. 
Therefore, combining all available data in the 0.2—10 keV range yields an RMS pulsed frac- 
tion of (36±2)%, which also remained constant over the entire course of the ~30 days of 
XRT observations. 

We also investigated the RMS pulsed intensity variations of SGR J1833— 0832 in time 
and energy using all available PCA data in the 2—4.5 keV, 4.5—8.5 keV and 8.5—14 keV 
intervals. These energy bounds are selected to accommodate enough counts in each interval 
and determine the RMS pulsed intensities with at least 5a significance. We find that the 



^The RMS pulsed fraction is defined as PFrms = (H(ElLi(Ri - ^avc)^ - ARf)^ ' / R^^,, , wliere N is 
the number of phase bins (N=f6), Ri is the rate in each phase bin, ARi is the associated uncertainty in the 
rate, and Rave is the average rate of the pulse profile. 
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Fig. 4. — {top) Plot of phase drifts for each segment of the PCA observations. The sohd hne 
is a second order polynomial model for the phase drift in time; {bottom) Fit residuals. 

pulsed count rates of the source remain fairly constant in all three energy intervals at around 
0.017, 0.025 and 0.006 counts/s/PCU, respectively. 



3.4. X-ray Spectral Properties 

We performed spectral analysis using only the XRT/PC observations, to avoid the high 
background effects of the WT observation. The data were rebinned to 20 counts per spectral 
bin to improve the fitting statistics. We fit the persistent source data with single component 
models and found that the first 12 PC mode observations (i.e., observations 000 through 
013 excluding 002 in Table [1]) are fully consistent with a constant source. Consequently, 
we summed all 12 observations into one single spectrum, which contained ~ 5000 counts 
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Fig. 5.— (left) Pulse profiles of SGR J1833-0832 in the 0.2-4 keV (filled circles), 4-6 keV 
(diamonds) and 6—10 keV (squares) energy bands; (right) normalized pulse profiles to the 
peak value of each energy interval. 



collected over 150 ks. We then rebinned the spectrum to 30 counts per bin and limited our 
spectral analysis to the 2 — 10 keV energy range, to account for the high absorption towards 
the source direction, and fit several absorbed (using phabs) single component models within 
XSPEC. The best fit spectral model parameters are listed in Table O Given the very high 
column density, the overall spectrum can be equally well fit by a steep power law (T = 3.4) 
or by a blackbody model with kT ~ 1 keV; no additional components are required. We 
note that the derived neutral hydrogen column densities (see Table |5]) are larger than the 
estimated Galactic column density of 1.7 x 10^^ cm~^ in the direction of SGR J1833— 0832. 
When we fixed the column density to this Galactic value, no single or double component 
model was able to account for the observed spectrum, indicating that the high column 
density could be intrinsic to the source or to its surroundings. It is important to also note 
that the origin of the high column density might be high density of molecular hydrogen in 
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the Galactic plane, on which the source is located. 

Table 5: Results of the spectral fits to the first 12 PC mode Swift/XBT observations. 

Model Column density T/kT 2-10 unabs.(abs) fiux x^ed 

(10^2 cm-2) -/(keV) (lO^^^ g^g ^^-2 g-i) ^^^f) 

Power law 16.4^};^ 3.4l°:2 17.3 (4.0) 1.03 (143) 

Black Body* 9.9+J]:^ l-H-aos 7.4 (3.8) 1.06 (143) 

Bremsstr. 13.71°;? 3.01°;^ 11.6 (3.9) 1.01 (143) 

*The equivalent black body radius for a source at 10 kpc is i?bb = 0.7 ± 0.1 km. 

We find evidence that the X-ray spectrum of SGR J1833— 0832 hardens as its fiux starts 
to decline in the last four PC mode observations (014 through 017 in Tabled]): an absorbed 
power law fit to these pointings yields indices in the interval from 2.6—3.2, and a black 
body fit yields temperatures in the interval from 1.2—1.3 keV (both at 90% confidence). 
The corresponding 2—10 keV unabsorbed (absorbed) power law and black body fiuxes are 
12.1x10"^^ (3.5x10"^^) erg cm"^ s"^ and 5.8x10"^^ (3.2x10"^^) erg cm"^ s"\ respectively. 

4. Discussion 
4.1. Burst emission 

At the onset of its active episode, SGR J1833— 08 32 emitted a single, relat ively faint. 



burst that was seen with Swift/BAT and INTEGRAL (JKuiper fc Hermsenll2010l ): the total 
event fiuence (15—350 keV) was 1.3x10"^ erg/cm^. Compared to the persistent emission 
fiux level (2—10 keV), the burst energetics are very low: it would take only 13 minutes for 
the persistent source to emit an equivalent amount of energy. It is, therefore, unlikely that 
the (detected) burst emission provided the energy budget for the persistent X-ray emission. 

The spectrum of the Swift/BAT burst is well described by a power law with expo- 
nential cutoff (Epeafc = 38 keV) or a black body with kT=10 keV. Its spectral and tempo- 
ral characteristics are very similar to those of known magnetars, in p articular, to those of 



SGR J0418+5729, which emitted only two rather dim bursts in one day (Ivan der Horst et al. 



20101 ). and IE 2259-1-586, which displayed a short episode of bursting for about 10 ks in 2002 



seen only in hard X-rays (JKaspi et al.ll2003l ). In addition to the Swift/BAT event, we iden- 
tified 4 weak X-ray bursts in our monitoring observations with the RXTE/PCA, which 
were not detected in the simultaneous Fermi/GBM data. It is interesting to note that the 
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RXTE/PCA bursts are very similar in duration and X-ray intensity to the five short and 
rather weak bursts, also detected in the RXTE/PCA passband, emitted by the young, rela- 
tively strongly mag netized (Brfjpo;e=4.9xl0^^ G) rotation powered pulsar PSR J1846— 0258 



flGavriil et al.ll2008[ ). 



4.2. Persistent emission properties 

The spin period, 7.5654091 s, and period derivative, 4.39x10^^^ s/s, of the source 
derived using the RX TE data are consistent (within 1.5cr) with the values obtained by 



Esposito et al. I (|2010[ ). Using these values, we estimate the dipole magnetic field strength of 
SGR J1833-0832 to be B = 1.8 xlO^^ G, which is at the lower end of the B-fields of SGRs 
and close to the median value of those of AXPs (Bmedian = 2.2x10^^ G). Note that the two 
recently discovered SGRs J0418 + 5729 and J0501 + 4516 possess the lowest SGR surface 
magnetic fields of < 3x10^^ G and 1.7x10^^ G, respectively. 

Figure |6]exhibits the evolution of the p ersistent X-ray e mission unabsorbed fiuxes (1 — 10 



keV) for five magnetars (data obtained from lRea et al.ll2009[ ) and SGR J1833— 0832, observed 



with Swift/XKT in the 2 — 10 keV range. Each source was observed immediately after an 
outburst, when the persistent emission level had significantly increased from its quiescent 
state. We notice that all sources had similar maxima at the onset of their activation, ranging 
between 10"^^ and 10"^° erg/cmVs. Three of the sources (SGRs J1627 - 41, J0501 + 
4516 and J1900 + 14) declined in a similar monotonic manner. However, SGR J1627 — 41 
settled at the lowest fiux level of all sources (~ 10~^^ erg/cm^/s), while SGR J1900 + 14 
reached a fiat plateau an order of magnitude higher and SGR J0501 + 4516 was still declining 
at the cessation of the XRT observations. The other two sources (IE 1547 — 5408 and 
CXOU J164710 — 455216) show similar trends (but different from the other three sources): 
they remain flat after activation for 10 — 20 days and decline thereafter with the same trend 
as the former three. 

We have plotted the SGR J1833— 0832 data as a band with an upper and lower bound 
on Figure |6l the upper flux limit is derived when we flt a single power law and the lower 
with a blackbody flt. Because of the excessive column density towards the source, there 
are very few photons below 2 keV and the spectral flts are not well constrained; the real 
unabsorbed flux of the source is within these two bounds, and seems to remain constant over 
the flrst ~20 days following the onset. Interestingly, we found an archival Chandra/ ACIS-S 
observation covering the fleld of SGR J1833— 0832 starting on 2009 February 13 10:50:20 
UT for a total of 8 ks. The source was not detected within our Chandra error circle ( §2.5p 
and we were able to derive a 2a upper limit (indicated by a dotted line on Figure E]) for 
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its quiescent unabsorbed flux of 3.4 x 10"^'^ erg/cm^ /s (lGogu§ Il2010bl ). This is almost 100 
times fainter than the current state of the source; assuming that SGR J1833— 0832 reaches 
this level again, it would take at least 2 years if it follows the same decline trend of the 
other five sources as indicated by the last XRT observations. The source quiescent level is 
close the lowest flux level of SGR 1627—41, the faintest magnetar thus far (Mereghetti et 
al. 2006, Esposito et al. 2009), and also to the flux of the faintest dim isolated neutron star 
RX J0420.0-5022 (Haberl et al. 2004) of 4.8 x lO-^^ erg/cmVs. 

Finally, pulsed radio emission has been detected from two magnet ars, XT EJ1810-197 
JCamilo et all l2006h and IE 1547.0-5408 fSGR J1550-5418: ICamilo et al.ll2007l ). The for- 
mer showed radio pulsations after an X-ray outburst from the source, but the latter was 
seen in radio while the source was quiescent in X -rays. Searches for pulsed radio ernission 
from other niagnet ars have been unsuccessful (e.g. iBurgay et al.ll2006t ICrawford et al.l 120071 : 
Hessels et al.l 120081 ). Our search for radio emission from SGR J1833— 0832 with the WSRT 
resulted in non-detections, with upper limits of 0.2 — 0.6 and 0.1 - - 0.3,mJy (dependin g on 
the assumed pulse widths) at 2280 and 1380 MHz, respectively. iBurgay et al.l (120101 ) re- 
ported an upper limit of 0.09 mJy at 1400 MHz on 2010 March 25 with the Parkes radio 
telescope. These upper limits are somewhat higher than some of the deepest radio searches 
for other magnetars, but more than an order of magnitude lower than the detections of 
XTEJ1810 - 197 and IE 1547.0 - 5408. 



4.3. Environment of SGR J1833-0832 



A visual inspection of the region surrounding SGR J1833— 0832 on the i^g— band images 
of the UKIRT GPS reveals no obvious source clustering in the vicinity of the source. How- 
ever, as the field is small and crowded, low mass open clusters may be difficult to discern 
from bright background sources. We have performed a color-color search to uncover any 
cluster members; we caution, however, that it is still likely that the number of "unrelated" 
sources may obscure the detection of any faint, low surface brightness cluster that may be 
present. 

There are several sources in the larger field which are plausible sites for the birth of 
SGR J1833-0832. One of these is the HH region IRAS 18310 - 0825, at RA = 18^33'^47^9, 

Dec = -08°23'52", which is offset by roughly 7' from SGR J1833-0832 another is [LPH96] 

023.1 62+0.023 {RA = 18''33"23^9, Dec = -08°40'32"; 10 arcmin offset. Ihockman. Pisano. fc Howard 
1996h . while a third is [KC97c] G023.5+00.0 at RA = 18^34'^20!ADec = ^8°21'23".7, 
which is also associated with a cluster known as BDS2003 118 ( jLockman Ill986[ ). Further, 
there are several dark clouds in the vicinity of SGR J1833— 0832, which apparently do not 
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yet host stellar populations. They have masses of 10^ — 10^ M©, and are likely locations for 
a next generation of star formation within the region. Finally, there is also a known pulsar 
(PSR J1830-08) roughly 37 away from SGR J1833-0832 with an adopted distance of 5.67 
kpc and a characteristic age of 147 kyr. It is possible that this source originates from the 
same star forming region associated with SGR J1833— 0832. 

While the dynamics of SGRs are as yet effectively unconstrained by observations (no 
SGRs have yet proper motion or radial velocity measurements) we can consider constraints 
by assuming that they i) have a short lifetime < 10^ years, and ii) have a kick velocity 
distribution broadly comparable to that seen for young pulsars. In practice this second 
assumption only improves the constraints we can make if the actual kicks are significantly 
smaller, and hence limit the transverse distance that SGR J1833— 0832 could have travelled 
since birth. A typical pulsar kick might be of order 200 km s~^, so that in 10*^ years the 
source may move ~ 2 pc on the sky. Hence the angular scale SGR J1833— 0832 could cover 
would be ^ ^ 400" /dkpc- For a distance of 6-12 kpc this corresponds to an angular scale of 
30 — 70", although allowing for a range of kick velocities might extend this range by a factor 
of 3 or more. In this case it would seem that none of the candidate clusters above are likely 
birthplaces. However, it is also possible (although not necessarily likely) that any progenitor 
may have "evaporated" from the cluster in which it formed with a velocity dispersion of a 
few km s~^. With a lifetime of a few Myr, this may result in a progenitor at a large distance 
from the cluster, making our constraints weaker. 
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Fig. 6. — Unabsorbed flux of SGR J1833— 0832 (filled circles) obtained by fitting power 
law (upper values) and blackbody (lower values) models to XRT data in PC mode, com- 
pared with other sources [SGR 1900+14 (plus), IE 1547.0-5408 (triangles), SGR 0501+4516 
(cros s), CXOU J164 710-455216 (diamonds) and SGR 1627-41 (squares); data obtained 



from lRea et al.ll2009| . Time is measured from the onset of burst active episodes. The hori- 
zontal dashed line indicates the Chandra 2a flux upper limit of SGR J1833— 0832 from 2009 
February 13. 



